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Glasma flux tubes vs minimum-bias jets in 2D angular correlations on rj and 
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Angular correlations measured in heavy ion collisions at the Relativistic Heavy Ion Collider 
(RHIC) include a same-side (SS) 2D peak which is strongly elongated on pseudorapidity rj in more- 
central Au-Au collisions. The elongated peak has been referred to as a "soft ridge." The SS peak 
is consistent with expected jet correlations in peripheral A-A and p-p collisions. A saturation-scale 
argument has been proposed to explain the origin of the elongated SS peak in terms of correlations 
from Glasma flux tubes interacting with radial flow. In this analysis we review the details of the 
proposed argument in comparison to perturbative QCD predictions of jet yields and correlations. 
We find that the proposal is inconsistent with several features of measured spectra and correlations. 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

Data from heavy ion collisions at the Relativistic 
Heavy Ion Collider (RHIC) have been interpreted within 
a hydro context to demonstrate formation of a thermal- 
ized, flowing partonic medium with small viscosity [H, 
However, differential spectrum and correlation analysis 
reveals a jet-like contribution whose variation with A-A 
centrality and collision energy seems to contradict hydro 
expectations In particular, a large-amplitude 2D 

peak at the origin in angular correlations on (77, </>) ex- 
pected from jet formation persists even in central Au-Au 
collisions, albeit the peak is elongated on pseudorapidity 
77 relative to a nominally symmetric jet cone 0, 

The elongated jet-like peak has recently been reinter- 
preted in terms of mechanisms other than parton scatter- 
ing and fragmentation, including "triangular flow" result- 
ing from initial-state transverse shape fluctuations (i.e. 
sextupole and other azimuth multipoles) which might 
in turn modulate radial expansion [6j and interaction of 
"Glasma flux tubes" with radial flow 0, H- 

In this article we confront the latter strategy, 
saturation-scale arguments based on Glasma flux tubes. 
We establish direct comparisons between saturation-scale 
predictions of particle production and measured yields. 
We examine a statistical argument for estimating the am- 
plitude of the same-side peak in terms of Glasma flux 
tubes coupled with radial flow. And we compare the i] 
dependence of theory predictions with measured 2D an- 
gular correlations as a test of saturation-scale theory. We 
conclude that the Glasma flux-tube description is falsified 
by spectrum and correlation data, whereas, perturbative 
QCD provides a quantitative description. 



A. 2D angular correlations 

Minimum-bias 2D angular correlations are constructed 
as autocorrelations on difference variables rj^ = rji — 772 
and 0A = - 02 |3|, y, 19|, ll^. No "trigger" parti- 
cle is used to define the event-wise angular origin. The 
only pt cut applied defines the lower limit of the momen- 
tum acceptance at pt = 0.15 GeV/c. So-called forward- 
backward correlation measurements are ID projections 
of 2D angular correlations onto 77A and so contain re- 
duced information. Dihadron azimuth correlations are 
ID projections onto 0a, with similar information loss. 

A per-particle correlation measure can be defined by 
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(1) 



where N is the mean multiplicity in angular acceptance 
(A77, A0) and histogram element rat is the ratio of event- 
wise sibling pairs to mixed pairs in 2D bin (a, b) within 
the pair angular acceptance. N{b)/Ar]A(j) = po{b) esti- 
mates the single-particle 2D angular density for impact 
parameter b averaged over the acceptance. It is related 
to the mixed-pair reference by po{b) w ^ Prej- 

The observed correlation structure in 2D histograms 
consists of a few components modeled by simple func- 
tional forms. Within the nominal STAR TPC acceptance 
(A?7, A0) = (2, 27r) the significant structure includes two 
?7A-independent multipoles (dipole and quadrupole), a 
same-side (SS) 2D peak, and a ID Gaussian on ?7a having 
negligible amplitude in more-central Au-Au collisions. 
The SS 2D peak and away-side (AS) azimuth dipole have 
been interpreted in terms of minimum-bias jets (minijets) 
with typical parton energy « 3 GeV hlHlJ. 



II. ANALYSIS METHOD 

We relate saturation-scale arguments and associated 
correlation predictions directly to measured particle pro- 
duction and minimum-bias 2D angular correlations from 
Au-Au collisions at y^J/vw = 200 GeV in the context of 
the Glauber model of A-A collision geometry. 



B. Glauber model of A-A geometry 

The Glauber Monte Carlo model relates collision pa- 
rameters Nhin (N-N binary collisions) and Npart (partici- 
pant nucleons) to collision geometry measured by impact 
parameter b and to a collision observable through cross 
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section a(b) determined via the minimum-bias distribu- 
tion on multiplicity rich within some angular acceptance. 

Centrality measure v = iNbm/Npart is a convenient 
measure for studying correlations relative to a Glauber 
linear superposition (GLS) reference in which A-A colli- 
sions are assumed to be linear superpositions of N-N col- 
lisions. The two-component model of A-A collisions de- 
scribed in Ref. '14] , when extrapolated from measured p-p 
collisions, is an example of a GLS reference. As a context 
for saturation-scale arguments note that R ^ N 



Ppart ~ Npart/R"^ ^ V and 7V( 
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III. SATURATION-SCALE ARGUMENTS 

Saturation-scale arguments applied to A-A collisions 
depend on longitudinal overlap of nucleon wave func- 
tions at small x [l^ and gluon transverse phase-space 
cell occupancy approaching unity. The result describes 
the distribution of gluons in a light-cone wave function 
produced by hard scattering off a large nucleus p^ . 



A. Saturation-scale gluon densities 

The problem addressed in Ref. (l6| is early-time ap- 
proach to thermalization, via multiple small-angle scat- 
tering in the transverse plane, of gluons released or 
freed from projectile nucleons near mid-rapidity. The 
saturation-scale treatment does not address large-angle 
(semihard) scattering of initial-state gluons. "At early 
times there is a negligible transfer of longitudinal mo- 
menta into transverse momenta" [l6| . In that descrip- 
tion gluons interact only with other gluons of the same 
rapidity during early-time thermalization. 

In Ref. [IGj b represents a transverse radius relative to 
the collision axis in central A-A collisions. The quark- 
gluon color factor is Cp = {N^ - l)/2Nc = 4/3. The 
strong-coupling constant is described by l/asiQ"^) ~ 
0.71n(QVA^), with A - 200 MeV [jj]. asCp is the 
coupling strength of valence (projectile) quarks to soft 
gluons. The projectile nucleon transverse density pN^b) 
integrates to A nucleons (for central collisions). The 
valence-quark density 3pAr(6) is the source for a field of 
soft gluons described in the low-density limit by 
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= 3pN{b)- 



ln(QVA^')(2) 



At high gluon densities saturation should occur when 
^Cpjcis gluons occupy each transverse phase-space cell. 
The saturation scale Q\ is then defined by 
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where c is the fraction of gluons released from projec- 
tile nucleons in central A-A collisions. Whereas in the 



low-density limit xG{x,Q'^) is slowly varying on Q^, in 
saturation xG{x,Qf) ^ l/as{Ql) ^ ln(Q2/A2). 



B. Saturation-scale hadron production 

We now consider the A-A centrality dependence of sat- 
uration following Ref. b represents the A-A im- 
pact parameter, and pN{b) — >■ Ppartib) is the participant- 
nucleon transverse density. Combining Eqs. ^ and ^ 
saturation sale Ql{x, b) is obtained by solving iteratively 
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TrasiQl)xG{x,Ql)ppart{b). (4) 



Because as{Q^) xG{x,QT) ^ 1 we expect Ql{x,b) ^ 

Ppart 

(6), and in Ref. ^ Ql ^ 2ppart/3 GeV^. But 
Ppart {b) K, v/2, where v = 2Nbin/Npart flol. So ~ ,y/3 
GeV^ , corresponding to the number of N-N binary colli- 
sions per participant pair. 

Integrating Eq. ([2]) over b^ (and y w ry) we obtain 



dNg{x,Ql) 
drj 



cxG{x,Ql)Npart{b) 



(5) 



with c^l. Assuming dNg/drj « (3/2) duch/dr] (LPHD) 
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Given l/as{Ql) - HQ'i/^'^), Ql ~ vji GeV^ and A 
0.2 GeV we obtain 
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where B is an 0(1) constant determined from data [l^l- 
For comparison, the eikonal Glauber two-component 
form is (first line) 



2 duch 
Npart dri 



ppp[l + xib)i,,-l)] (8) 
Snn + vHAAib), 



with Pp-n_~ Snn + Hnn, the peripheral limit of A-A col- 
lisions H HI ll!]. Parameter x{b) « HAA{b)/{SNN + 
Hnn) (second line) actually has a strong centrality de- 
pendence, in contrast to the fixed x value in Ref. [14l|. 



C. Comparison with hadron production data 

We can compare the saturation-scale prediction with 
measured particle production and angular correlation 
trends. In Fig. [1] (left panel) the jet-correlated pair den- 
sity (same-side 2D peak integral) is averaged over a lim- 
ited angular acceptance (solid curve) The dashed 
curve shows the result for An angular acceptance. By 
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FIG. 1: Left: Centrality dependence of the same-side 2D peak 
volume (solid curve) measured by pair ratio j^{b) within ac- 
ceptance 77 £ [—1, 1] in combination with single-particle an- 
gular density po{b). Right: Particle production trends vs cen- 
trality f from spectrum data (points), from minijet angu- 
lar correlations (bold solid curve), from the eikonal Glauber 
model of Ref. with x = 0.09 (dash-dotted line) and from 
a saturation-scale argument (bold dashed curve). 

invoking pQCD jet number Nj{b) per A-A collision ap- 
pearing within the angular acceptance the correspond- 
ing jet fragment density can be inferred, represented by 
quantity vHAA{b) [Tsj . 

Two-component particle yield model (2/Npart)po{b) — 
Ppp[l + x{v — 1)], representing the eikonal model in 
Ref. :14,] with x = 0.09, is shown as the dash-dotted 
line in Fig. [T] (right panel). Sohd curve Snn + vHAA{b) 
includes HaaW) inferred from two-particle jet correla- 
tions in Ref. '20'] which increases by factor 4.5 from pe- 
ripheral to central Au-Au collisions, factor 1.5 represent- 
ing increase of the dijet cross section and factor 3 rep- 
resenting increase of the mean jet fragment multiplicity 
arising from fragmentation function modification 13, 20]. 
The effective x{h) value (slope on v) thus increases from 
X « 0.02 inferred directly from p-p spectra [l^ (trend 
labeled GLS) to x « 0.1 consistent with more-central 
Au-Au yields. 

The bold dashed curve is saturation-scale trend 
(2.5/2.25) 0.82 ln(8i^) from Eq. (19) of Ref. Q scaled 
from 130 to 200 GeV and using the equivalent of Q^/A^ 
determined in the present analysis. Prefactor B = 0.82 
established in Ref. [l3j matches 130 GeV Au-Au data 
in more-central collisions. If the prediction is scaled to 
match more-central yield data (where saturation would 
be more likely) it fails dramatically for more-peripheral 
collisions. The prediction is concave downward, whereas 
the spectrum data (points [5]) are significantly concave 
upward. 

The two-component models in Fig. [T] (right panel) as- 
sume soft component Snn ~ 2.4 independent of cen- 
trality, any yield centrality dependence arising from par- 
ton scattering and fragmentation, whereas saturation- 
scale arguments predict a soft-component yield varying 
strongly with centrality as ln(8i/) with no significant par- 
ton scattering contribution. The soft (non-jet) compo- 
nent can be divided conceptually into a non-saturation 
part dominating p-p particle production (conventionally 



described by string fragmentation and/or soft Pomeron 
exchange) and a conjectured saturation part depending 
strongly on centrality but with no prediction of absolute 
magnitude. Because addition of fixed soft component 
Snn "^o measured hard component vHAAib) derived in- 
dependently from jet correlation data matches the spec- 
trum yields (points) we conclude that the flux-tube con- 
tribution is negligible even in central Au-Au collisions. 

IV. 2D ANGULAR CORRELATIONS 

Minimum-bias (pt-integral) 2D angular correlations 
provide an essential reference system for any correlation 
analysis in which pt cuts are imposed to study jet phe- 
nomenology (e.g. trigger-associated analysis) or azimuth 
quadrupole systematics. 2D correlation data are accu- 
rately described by a simple fit model including three 
principal model components: (a) a same-side 2D peak 
at the origin well approximated by a 2D Gaussian for 
all minimum-bias data, (b) an away-side ridge well ap- 
proximated by AS azimuth dipole [1 — cos((/)a)]/2 for 
all minimum-bias data and uniform to a few percent 
on TjA (having negligible curvature), and (c) an azimuth 
quadrupole cos(20a) also uniform on t^a to a few percent 
over the full angular acceptance, (a) and (b) together 
have been interpreted as minimum-bias jets or "mini- 
jets" [1^. (c) is conventionally interpreted to represent 
elliptic flow, a hydrodynamic phenomenon. The 2D data 
are described to a few percent of the major-feature am- 
plitudes over the entire angular acceptance. Fit residuals 
typically contain no significant nonstatistical structure. 




FIG. 2: Left: Amplitude of a same-side 2D Gaussian fitted 
to minimum-bias 2D angular correlation data from 200 GeV 
Au-Au collisions 4]. Right: Fitted peak widths for the same- 
side 2D Gaussian. GLS indicates a Glauber linear superposi- 
tion reference extrapolated from measured p-p collisions [l^ . 

Figure [5] summarizes preliminary fitted SS 2D peak pa- 
rameters vs centrality measure v within nominal STAR 
TPC angular acceptance (A77, A0) = (2,27r) 3- ^20 in 
the left panel is the fitted amplitude of the SS 2D Gaus- 
sian function. Its two r.m.s. peak widths are reported in 
the right panel. There is smooth variation with central- 
ity, but a "sharp transition" in SS 2D peak properties 
occurs at a specific point on centrality — v ~ 2.5. The 
corresponding quadrupole data are reported in Ref. [2l| . 




FIG. 3: (Color online) Angular correlation histograms for four 
centralities of 200 GeV Au-Au collisions based on fit param- 
eters from Ref. 3|. The centralities starting from upper left 
are given hy v = 1.25 (- N-N ), 2.5, 4.5 and 6 (6 = 0). 



Figure [3] shows examples of 2D angular correlations 
from four centralities of 200 GeV Au-Au collisions within 
acceptance ry G [—1,1] based on the fit parameters of 
Ref. Q. The centralities correspond to v = 1.25 (~ N-N 
collisions), 2.5, 4.5 and 6 (6 = 0). The histograms are 
plotted within the STAR TPC angular acceptance usu- 
ally adopted for 2D correlation analysis. 

Figure [3] (upper right) with v = 2.5 (the transition 
point) includes a SS 2D peak with 1:1 aspect ratio (both 
peak widths 0.8). Because the plotting format is 
27r : 4 ^ 3 : 2 the symmetric peak seems to be elon- 
gated on 77. Apparent peak elongation should be checked 
against the actual data (model fit) properties and plot- 
ting format. The SS peak seems to be superposed on 
a distinct "ridge" structure due to apparent peak elon- 
gation and the underlying azimuth quadrupole contribu- 
tion. Those data can be compared with Fig. 1 of Ref. [1^ . 

Arguments in Refs. 0, ll| pertain to the SS 2D peak 
at the angular origin. Figure S] (left panel) illustrates 
nominal jet angular correlations in more-central Au-Au 
collisions (v — 4.5, lower-left panel in Fig. [3]) revealed 
when non-jet model components (mainly the azimuth 
quadrupole) are subtracted from 2D data histograms. 
The remaining elements are the SS 2D peak and the AS 
dipole. Figure [4] (right panel) shows the isolated SS 2D 
peak (dipole subtracted) extrapolated to larger 77 accep- 
tance 77 S [—1.5, 1.5], which is actually accessible within 
the STAR TPC albeit with increased systematic uncer- 
tainty. The SS peak for minimum-bias angular correla- 
tions is consistent with a single 2D Gaussian in more- 
central Au-Au collisions. SS peak structure retains a 
large curvature on 77A in all cases. 



FIG. 4: (Color online) Left: Angular correlations from 200 
GeV Au-Au collisions with centrality v = 4.5 extrapolated 
to G [—1.5, 1.5] with nonjet fit components (quadrupole, ID 
Gaussian on t^a) subtracted to reveal nominal jet correlations. 
Right: The previous histogram with the away-side dipole term 
subtracted to isolate the same-side 2D jet peak. 



V. GLASMA FLUX TUBES AND THE SS PEAK 

Variation with A-A centrality of the amplitude and 
azimuth width of the SS 2D peak is described in terms 
of Glasma flux tubes, which are said to be "pushed out" 
by radial flow to form a SS ID peak narrow on azimuth 
and uniform on 77A . A statistical argument is introduced 
to express peak amplitude scaling in terms of saturation 
scale Q\. A blast-wave model is invoked to determine the 
effects of radial flow on peak azimuth width and possibly 
amplitude variation. 



A. Statistics of a two-tiered hierarchy 

The treatment of correlations from Glasma flux tubes 
in Ref. 0] invokes a cluster or clan model [11] . Such mod- 
els apply as well to longitudinal string fragmentation |23| 
and transverse-scattered parton fragmentation 13, 201. 
Number fluctuations within angular acceptance (A?7, A(/)) 
correspond to a running integral of angular correlations 
up to those limits [§] . H. denotes a per-pair integral fluc- 
tuation measure comparable to correlation measure r — 1 
deflned in Ref. [3|]- A corresponding per-particle integral 

measure is NTZ = '^"^^^ — Aaj^{Ar], A(j)), where the last 
quantity is a variance difference — defined as a running in- 
tegral (on angle bin size or scale) of angular correlations 
in the per-particle form A/9(?7a, 0A)/^/pre7 0, US- The 
integral equation can be inverted to obtain differential 
number angular correlations within the angular accep- 
tance from fluctuation scale dependence @ . Inversion of 
mean-pt fluctuation scale dependence led to identification 
of minijet pt (as opposed to number) angular correlations 
as the source of those fluctuations [Ifl, [23 . 

If gluon/hadron number N results from event-wise pro- 
duction of K sources, each emitting n particles (both are 
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random variables fluctuating about mean values) then 



N 



N 



+ A- 



'K 



K 



K 



.(9) 



illustrating the additivity of per-particle measures. That 
expression could describe fragmentation of flux tubes, 
prehadrons from a Lund string or large-angle scattered 
gluons. The first term is a variance difference represent- 
ing the integral of intra-source correlations — the SS 2D 
peak. The second term is a covariance representing the 
integral of source-source correlations — the AS ridge. A 
measures the fraction of sources (e.g. jets) with a single 
partner in the 77 acceptance, anticipating a back-to-back 
jet interpretation. The third term represents event-wise 
source-number fluctuations due to A-A centrality fluctu- 
ations and centrality fluctuations of N-N collisions within 
A-A collisions. The fourth term represents a contribution 
due to the two-tiered hierarchy mechanism alone, present 
even when sources and fragments are Poisson distributed. 
The third and fourth terms represent correlation contri- 
butions uniform on angle, not resolved peak structure. 

In Ref. Jl the first and third terms of Eq. ([9]) are as- 
sumed to be negligible (Poisson processes), and the co- 
variance term is not acknowledged. The fourth term n 
is attributed to the SS 2D peak amplitude. NTZ — n or 

The mean num- 



7^ - l/K, with K = QlR^ 
ber of fiux tubes (sources) is proportional to the number 
of participant nucleons. Since N ^ [l/as((3^)] Q^i?^, 
idN/dr])n - n - l/as(Q?) \n{8iy) is taken as a mea- 
sure of the mean per-particle 2D peak amplitude. That 
is also the trend claimed for the per-participant-pair par- 
ticle yield {2 / Npart) duch / drj in Ref. |!l4], as in Eq. 
Because the fourth term of Eq. © describes a uniform 
background contribution to angular correlations it is not 
relevant to the SS 2D peak. 

In a jet interpretation of Eq. ^ the first term on the 
RHS measures intra-jet correlations — the SS 2D peak 
integral (not the amplitude), which is the actual sub- 
ject of Ref. 0. The total number of correlated pairs is 
K{a^ — fi) — )■ Nj{b) n^^ -{h), a product of the mean event- 
wise jet number and the number of fragment pairs per 
jet. A pQCD calculation of Nj{b) combined with mea- 
sured angular correlations from Ref. [3| quantitatively 
describes single-particle spectrum hard-component yields 
attributed to jets, as in Fig. [1] (right panel) [20|. FigureH] 
(left panel) shows two of the four elements contributing 
to NTZ in Eq. (0). The integral of the SS 2D peak is 
represented by the first term in the RHS. The AS dipole 
peak integral is represented by the second term. 
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B. Glasma flux tube centrality trend 

In Ref. [7] the SS 2D peak amplitude for A-A collisions 
is defined by 



y/Pref 



na^ ^[Qs{^Jsnn, b)] F{(j)^; ^snn, b),{10) 



where -F'((/>a) is a unit-normal peaked function centered 
at zero azimuth difference, its width ^ 0.65 varying only 
slightly with energy and centrality. The amplitude of 
F{(f>) cx can then vary by at most 20% given the 

slowly- varying trend of cr^^ in Fig. 2 of Ref. 7]. How- 
ever, F(6^) is attributed to the dashed curve in Fig. 1 
of Ref. [7| which varies by more than a factor 3. 

In Ref. the strong-coupling constant is given by 



l/a.iQs) « 0.71n(g2/o.04GeV2), 



(11) 



with Q2 ~ GeV^ from Ref. [Tl|. Thus, the SS 2D 
peak amplitude is predicted to vary as ln(8i^), just as 
for saturation-scale particle production in Eq. ([7]). The 
peak amplitude prediction for 200 GeV Au-Au in Fig. 1 
(upper panel, solid curve) of Ref. [Tj] varies between 0.05 
and 0.75, increasing by a factor 15 over 1 < ly < 6. But 
ln(8 X 6)/ln(8) = 1.9. The solid curve from Ref. [7] then 
seems to be incompatible with Eq. ([TU|. 

The aj^ factor in Eq. IjlOp is derived from an argument 
based on fluctuation measurements which would relate 
to the volume of a bounded 2D correlation peak, not the 
amplitude of a ID peak on azimuth unbounded on rj. 



C. Glasma flux tube correlation description 

Figure [5] provides a direct comparison between the 
full flux-tube correlation prediction and measured data. 
In Figure [5] (upper panels) the flux-tube prediction is 
compared to peripheral Au-Au « p-p angular correla- 
tions {ly = 1.25 in Fig. [3]) extrapolated to 77 acceptance 
77 € [—2, 2]. The ID Gaussian component on rjA has been 
removed, leaving SS and AS jet contributions. Evaluat- 
ing the flux-tube ridge structure as a possible component 
of p-p correlations we conclude that the 2D data impose a 
small upper limit on ridge structure consistent with zero. 

In Figure [5] (lower panels) the SS 2D peak in data from 
central Au-Au collisions (1/ = 6, the quadrupole is neg- 
ligible) is substantially elongated on 77A. However, the 
2D data (left panel) do not require a significant fiux- 
tube ridge structure (right panel). Although some trig- 
gered 2D angular correlations seem to suggest a separate 
"ridge" extending beyond a localized and symmetric 2D 
"jet" peak [2^ the selected particle pairs are a subset of 
minimum-bias data which do not so indicate. If pt cuts 
are imposed significant deviations from an ideal Gaus- 
sian shape should be expected due to local charge, fla- 
vor, baryon-number and momentum conservation during 
parton fragmentation. For instance, SS peak shapes for 
like-sign and unlike-sign pairs may be quite different j27| . 



VI. DISCUSSION 

The principal goal of Refs. 0, IB| is reinterpretation of 
minijet phenomenology (e.g. systematics of SS 2D peak 
in angular correlations) in terms of a non-jet mechanism 
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FIG. 5: (Color online) Comparison of fits to 2D data his- 
tograms (left panels) with Glasma flux tube predictions (right 
panels) for Au-Au collisions with v = 1.25 (« N-N collisions, 
upper panels) and v = % (h = lower panels). 



based on Glasma flux tubes. We review major elements 
of the Glasma flux-tube argument in comparison with a 
pQCD description of minijet manifestations in spectra 
and correlations. 



A. Ridge terminology and r\ dependence 

Imposition of the term "ridge" on SS 2D correlation 
structure has caused significant confusion. The SS corre- 
lation structure is a single 2D peak under almost all con- 
ditions, although some -pt cut conditions produce signifi- 
cant deviations from an ideal Gaussian shape. The com- 
bined terms "soft" and "ridge" seem to distance SS corre- 
lation structure from a nominally-ideal jet phenomenon 
manifested by high-pt (hard) particles falling within a 
symmetric jet cone. 

The implied dichotomy of a hard (jet) peak and soft 
(non-jet) ridge is questionable. Most parton fragments 
in pQCD jets appear below 2 GeV/c O, [13, just where 
thermalization and hydro flows are commonly believed 
to dominate RHIC collisions 28*]. Such "soft" fragmen- 
tation is susceptible to strong angular deformation de- 
pending on the collision system. The SS jet peak in p-p 
collisions is strongly elongated (2:1) in the azimuth di- 
rection [2^ in contrast to comparable elongation (3:1) in 
the "T] direction in central Au-Au collisions [3, 0, [2^ . 

As noted in Sec. IV C[ accurate measurement and de- 
scription of correlation r\ dependence is essential for cor- 
rect physical interpretation. 77 dependence permits sep- 
aration of 2D correlations into ?7-independent (negligible 



curvature on ?7a) and 7y-dependent (strong curvature on 
?7a) structure within the STAR TPC acceptance. 



B. Glasma flux tube correlation mechanism 

According to Refs. 0, [1| early-stage rapidity corre- 
lations combined with late-stage radial flow produce a 
SS correlation structure elongated on 77A. Early-stage 
correlations are attributed to fragmentation of Glasma 
flux tubes to separate gluons. Late-stage radial flow 
results from thermalization of the same flux-tube frag- 
ments. Both the proposed correlation mechanism and 
persistence of any such correlations in the measured final 
state are questionable for four reasons: 

a. Absence of conjectured initial correlations 
The concept of (color) correlation length seems to be 
confused with isolated flux tubes acting as distinct par- 
ticle sources. In the case of isolated "string" fragmenta- 
tion in ip-p collisions (flux tube diameter ~ hadron size) 
described by the Lund model we observe angular cor- 
relation structure narrow on ?7a and nearly uniform on 
(fiA m, [13] • That (soft component) correlation structure 
falls to zero amplitude by mid-centrality in Au-Au colli- 
sions The concept of correlated gluon emission from 
individual Glasma flux tubes is a theoretical conjecture 
unconflrmed by experimental evidence. 

b. Time sequence: boost vs emission Accord- 
ing to Ref. particles (gluons) resulting from flux tube 
fragmentation thermalize and produce radial pressure 
gradients which drive radial flow. The radial flow then 
boosts the correlated particles to produce narrow struc- 
tures on azimuth. We might expect a large correlation 
effect when particles are emitted from a common source 
already boosted, in which case isotropic emission in the 
boost frame should appear as a directed jet of particles 
in the lab frame. A notable example is QCD jets. 

If particles are emitted isotropically from a source sta- 
tionary in the lab and then boosted radially, some pair 
opening angles will increase while others decrease. We 
should expect a minimal correlation effect. Thus, the 
time sequence is critical. For the sequence described in 
Ref. [7j , even if flux-tube fragments were initially strongly 
correlated in conflguration space the thermalization pro- 
cess (which would also tend to reduce initial-state cor- 
relations) followed by late radial boost should result in 
negligible momentum angular correlations. 

c. Absence of radial flow Radial flow is inferred 
from pt spectra by blast-wave flts within limited pt in- 
tervals below 2 GeV/c, based on the assumption that 
the entire spectrum in that pt interval corresponds to 
emission from a thcrmalized bulk medium. However, 
differential analysis of 200 GeV Au-Au pt spectra for 
two hadron species reveals that spectrum evolution is 
dominated by the spectrum hard component, with most- 
probable Pt ~ I GeV/c Blast- wave parameters in- 
ferred from spectrum flts correspond closely to fragment 
distributions predicted quantitatively via pQCD [l3l[28|. 
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d. Central-limit attenuation According to the 
argument in Ref. Q almost all final-state hadrons come 
from flux tubes. The number of flux tubes in central 
Au-Au coUisions is Q'j.R'^ « Npart ~ 200. Flux tubes 
at different radii should suffer different radial boosts and 
thus correspond to different peak widths. All such con- 
tributions would then be superposed event-wise and ran- 
domly distributed on the 27r azimuth acceptance. Super- 
position should lead to what can be termed central-limit 
attenuation: the many contributions sum to a near con- 
tinuum on (j) with small-amplitude fluctuations. 



C. Glasma flux tubes compared to QCD jets 

According to Ref. [7] Glasma flux tubes effectively re- 
place minimum-bias jets as the fundamental QCD pro- 
cess in nuclear collisions. The statement "Flux tubes 
arise naturally in... high energy collisions" does not differ- 
entiate between Lund strings well-studied in elementary 
hadronic collisions and conjectured Glasma flux tubes. 

Flux tubes and large-angle scattered partons do share 
some formal similarities. Both objects may fragment to 
daughter gluons which then convert to hadrons according 
to LPHD. Both might result in SS angular correlations 
among detected hadrons. Both are nominally described 
by QCD, one in the perturbative limit, the other in the 
continuum limit. Both might involve transport from lon- 
gitudinal to transverse momentum space. But the two 
possess very different systematics which disfavor the for- 
mer in comparisons with data. Minimum-bias jets con- 
form to all aspects of pQCD in spectrum and correlation 
measurements except rj broadening, whereas Glasma flux 
tubes conform to almost no measurement features except 
the rj broadening which motivated their introduction. 

a. Glasma flux tubes In this one-component 
model the (soft) component arises from longitudinal frag- 
mentation of projectile nucleons (via flux tubes). Large- 
angle parton scattering and fragmentation do not play 
a significant role. Flux tubes are transverse phase-space 
cells in an 77-independent theory. The single parameter 
of the theory is the saturation scale oc v. l/as{Ql) 
said to account for all y/s and centrality dependence. 

The multiplicity of flux tubes is ^ ^ Q^R^ oc 
Npart- The mean fragment multiplicity per flux tube 
is n ~ 1/asiQl) ^ HQl/^'^) ^ HM- The to- 
tal gluon/hadron multiplicity from all flux tubes in an 
A-A collision is Kn - Q'^^R\\n{Ql/A'^) - NpartMM- 
The per-particle SS correlation amplitude is predicted 
by (drich / dr])TZ — n ^ InlSv). There is no AS peak 
on azimuth in the theory — no tube-tube correlations, no 
momentum transfer between flux tubes. 

The flux-tube system relies on conjectured radial flow 
to translate longitudinal momentum to transverse phase 
space (via thermalization and radial pressure gradients) . 
Inference of radial flow from spectra relies on interpreting 
parton fragmentation as a hydro manifestation 0, [2^ . 



b. Minimum-bias jets In the two-component 
model [1^ [1^ participant nucleons fragment longitudi- 
nally to fbced hadron multiplicity (Lund model) indepen- 
dent of A-A centrality (soft component). The contribu- 
tion from minimum-bias jets (minijets) increases at least 
as fast as Ntin, rising to 1/3 of the total rich in central 
Au-Au collisions at 200 GeV [20] (hard component). 

The multiplicity of minimum-bias (~ 3 GeV) jets in 
acceptance A77 is Nj — Ni,inArjf{b), where N-N jet fre- 
quency f{b) is slowly varying with centrality [itI, |2^. 
The fragment multiplicity periet in p-p collisions is 
nch.j ^ ylnax ~ [log(Q/m^)]2 [l3], where Q « 6 GeV is 
the dijet energy, and increases to three times that num- 
ber in central Au-Au collisions, consistent with strong 
modiflcation of fragmentation functions without signifl- 
cant net energy loss from jets p^. Total multiplicity rich 
in A-A collisions is the sum of soft and hard components. 

For angular correlations the per-particle SS 2D peak 
volume is measured by NTZ — Nj n^i^ j/rich [IQi]- Angu- 
lar correlations also include a prominent AS ID peak on 
azimuth (AS ridge) whose centrality variation (cx Nun) 
closely matches the SS 2D peak and is very different 
from that expected from global momentum conservation 
(oc Npart)- AS ridge characteristics require large-angle 
scattering of parton pairs, with each parton localized on 
X and with a broad distribution on difference x-\ — X2 , not 
longitudinally co- moving as assumed in Ref. [16|. 

A specific energy dependence is expected for jet phe- 
nomena cx log{y/SNN GeV} (inferred from V2\2D} 
data obtained from fits to 2D angular correlations 21J). 
According to that trend the 62 GeV SS 2D peak am- 
plitude should be a fraction 0.57 of the 200 GeV peak 
amplitude for Au-Au collisions, and that is observed [ij. 

c. Comparison summary The minijct descrip- 
tion has good predictive power for hadron spectrum and 
correlation structure. The spectrum hard component is 
predicted quantitatively in terms of minimum-bias frag- 
mentation by folding a pQCD dijet spectrum with mea- 
sured fragmentation functions (modified in more-central 
A-A collisions). [H, [2^. The spectrum soft component 
may be a universal feature of all collision systems [28| . 

Event-wise reconstruction of jets in p-p collisions has 
been achieved down to 3-4 GeV by UAl [3l| as well as 
STAR [s^l- Minijet systematics in Au-Au collisions in- 
ferred from particle production, spectra and minimum- 
bias correlations agree quantitatively with results from 
jet reconstruction in terms of (a) no jets lost to ther- 
malization, (b) particle-number suppression at larger pt 
counterbalanced by larger enhancement at smaller pt cor- 
responding to (c) nearly conserved total jet energy. En- 
ergy is rearranged within jets and angular correlations 
are modified by SS peak elongation on rj, in contrast to 
elongation on cj) for p-p collisions. 

The description based on Glasma flux tubes fails to 
provide a comparably accurate and comprehensive de- 
scription. Theories which ignore 2D correlation structure 
on 77 (e.g. curvatures), which deal only with projections 
onto azimuth, and which do not account for differential 
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Pt spectrum structure and its evolution with centrality 
and collision energy are not adequate to describe the en- 
semble of experimental results now available from RHIC. 



D. Flux tubes, causality and cosmology 

It is argued in Ref. 0] that "ridge" correlations ex- 
tending over large rj intervals require a source established 
at early times. "Correlations over several rapidity units 
can only originate at the earliest stages of an ion col- 
lision when the first partons are produced. ...causality 
limits such effects [hydro modifications] to a horizon of 
roughly from one to two rapidity units." The argument 
can be summarized as (a) "Flux tubes arise naturally 
in... high energy collisions." (b) Glasma flux tubes form 
very early in the collision (< 1 fm/c) (c) long-range corre- 
lations require an early-time source (causality argument), 
therefore (d) long-range correlations must be produced 
by Glasma flux tubes. 

Regarding (a,b) - Two types of flux tubes can be dis- 
tinguished. The flux tubes described in [13] of Ref. Q 
(e.g. Lund model) are indeed expected in QCD but have a 
nonperturbative transverse size measured by 1 / Aqcd ~ 
1 fm, whereas conjectured Glasma flux tubes are said to 
have a perturbative transverse size measured by 1/Qs ~ 
0.1 — 0.2 fm. The latter are not an inevitable consequence 
of QCD as implied by Ref. 0. 

Regarding (c,d) - 77 is a measure of polar angle. Frag- 
mentation extending over large time intervals can pro- 
duce signiflcant yields at polar angles far from tt/2 and 
hence at large rj, especially with the large increase of 
small-pt jet fragments in more-central Au-Au collisions. 

An alternative broadening mechanism based on con- 
ventional pQCD is color connection of the scattered par- 
ton to a projectile nucleon following hard Pomeron (color 
singlet) exchange (see Sec. XIII-C,D of Ref. [13]). That 
mechanism naturally produces jet elongation in one di- 
rection on 77 which would be symmetrized in an ensemble 
average. It is also an "early time" phenomenon capable 
of producing structure over large rapidity intervals. 



E. Pt dependence of correlated pairs 

In Fig. 5 (lower panel) of Ref. [1] correlation data 
are shown for pt cuts which correspond to a running in- 
tegral from below of the SS 2D peak on pt- The data 
as presented give the impression that the most proba- 
ble particle momentum in the SS 2D peak is 0.2 GeV/c, 
apparently supporting the imposed "soft ridge" termi- 
nology. However, the underlying spectrum of correlated 
particles corresponds to the negative derivative of the 
plotted data. The true most-probable momentum for 
those data is 1 GeVJc, consistent with spectrum hard 
components for p-p [l9| and more-peripheral Au-Au col- 
lisions [Hj . Spectrum hard components arc in turn quan- 
titatively described as jet fragment distributions by the 



folding of a pQCD dijet spectrum with measured frag- 
mentation functions [ij]. The term "soft" is misleading. 

The pttrend of the SS peak azimuth width is described 
in Ref. 8] in terms of radial flow: larger pair pt corre- 
sponds to smaller opening angle because a fluid cell has 
larger radial boost. But the same trend is expected from 
pQCD jets, with "fluid cell" translated to location in the 
fragmentation cascade descending from a parent parton. 
That is one of several examples where pQCD jet phe- 
nomenology is recast within a hydro scenario. Another 
is interpretation of low-pf jet fragment distributions in 
terms of radial flow [isj . Distinctions made between a 
"thermal bulk" at smaller pt and jet fragments at larger 
Pt are unjustified, since according to pQCD (MLLA |33l ]) 
most par ton frag ments appear below 2 GeV/c for any 
jet [sl-lTsl. [itI . [28| , also implying that most jet-correlated 
hadron pairs appear at larger angular separations. 



F. Spectrum and correlation references 

For differential spectrum and correlation analysis of 
A-A collisions it is essential to define a Glauber linear 
superposition (GLS) reference based on differential mea- 
surements in elementary collisions. The GLS reference 
for two-component spectrum analysis is [1, [l^] 



iV, 



part 



= SNNiVt) + i^HNNivt) (12) 



reflecting binary-collision scaling of fixed hard compo- 
nent Hnn. In Ref. [4| the dashed correlation reference 
curve in Fig. 3 (left panel) represents the GLS reference 



0M5NMn 



0.045 



-[SNN + yHAA{h)]/Ppp l + x{b){v~l) 



(13) 



extrapolated from the SS 2D peak in 200 GeV p-p colli- 
sions, where x{h) increases linearly from 0.02 in p-p colli- 
sions to X ^ 0.1 in central Au-Au collisions, x — 0.02 « 
H]sin/{Snn + Hnn) in acceptance A?/ = 2 is derived 
from a two-component analysis of p-p pt spectra p^ . 
X ^ 0.1 matches spectrum yields in more-central Au-Au. 

In Ref. 8] similar dashed curves in Fig. 3 are described 
as "wounded nucleons + fiow." "Flow" should play no 
role in extrapolations from p-p collisions. The preferred 
(GLS) reference is based on known binary- collision scal- 
ing of measured jet correlations [l^i and spectrum 
trends [l9l] in p-p collisions. 



VII. SUMMARY 

A saturation-scale argument has been proposed to ex- 
plain the origin of the ry-elongated same-side (SS) 2D 
peak ("soft ridge") in pt-integral (minimum-bias) angu- 
lar correlations in terms of Glasma flux tubes plus radial 
flow. The SS peak is consistent with expected number 
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and pt angular correlations from jets in more-peripheral 
A-A and p-p collisions. The ?7-elongation phenomenon 
has been observed in more-central Au-Au collisions. 

In the present analysis we have reviewed saturation- 
scale predictions for particle production and find that 
the prediction oc log(8i^) [v = 2Npart/Nbin measures the 
number of binary N-N collisions per participant pair) dis- 
agrees with measured hadron production, whereas a two- 
component pQCD analysis of jet angular correlations de- 
scribes the hadron production data quantitatively. 

We have examined a statistical argument relating the 
saturation scale to the amplitude of the SS 2D peak and 
find that two contributions to number fluctuations from 
angular correlation structure have been confused. The 
contribution which might correspond directly to the vol- 
ume of the SS 2D peak (intra-source correlations) is as- 
sumed to be negligible, whereas a contribution which 
corresponds to a uniform correlation background is at- 
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